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ABSTRACT 

The dominant component of the (100 MeV - 50 GeV) GRB emission detected by LAT starts with 
a delay relative to the prompt soft (sub-MeV) gamma-rays and lasts long after the soft component 
fades. This has lead to the intriguing suggestion that this high energy emission is generated via 
synchrotron emission of relativistic electrons accelerated by the external shock. Moreover, the limits 
on the MeV afterglow emission lead to the suggestion that, at least in bright GeV bursts the field is not 
amplified beyond compression in the shock. We show here that considerations of confinement (within 
the decelerating shock), efficiency and cooling of the emitting electrons constrain, within this model, 
the magnetic fields that arise in both the upstream (circum burst) and downstream (ejecta) regions, 
allowing us to obtain a direct handle on their values. The well known limit on the maximal synchrotron 
emission, when combined with the blast wave evolution, implies that late photons (arriving more than 
^ 100 s after the burst) with energies higher than ~10GeV do not arise naturally from external shock 
synchrotron and almost certainly have a different origin. Finally, even a modest seed flux (a few mJy) 
at IR-optical would quench, via Inverse Compton cooling, the GeV emission unless the magnetic fleld 
is significantly amplified behind the shock. An observation of a burst with simultaneous IR-optical 
and GeV emission will rule out this model. 
Subject headings: gamma rays: bursts 



1. INTRODUCTION 

The recent observations of the Large Area Telescope 
(LAT) on board of Fermi of GeV emission (lOOMeV 
- 50GeV) from GRBs revealed an interesting pattern. 
The GeV emission is del ayed relative to the o nset of the 
prompt MeV emission (Abdo et al. 2009a). It shows 
a constant po wer-law decay long affer the prompt emis 



sion dies out (l Abdo et al. |2009bt [Ghisellini et al |2010[ ). 

While surprising at first, one may recall that a "precur- 
sor" of these observations was made already by EGRET 
that detected an 18 G eV photon 90 minu tes after the 



burst in GRB 940217 Hurley et al. (19941 an d a rising 



late GeV spe ctral component m GRB 941017 ( Gonzalez 



et al. 



20031). This pattern suggests that the bulk of 
the GeV eniission arises from an external shock after- 
glow pCumar fc Barniol Duran||2009a|b Ghisellini et al 



2010[ r While a detectable high ener gy external shock 
emission was expected for a long time dMeszaros & Rees 



[1994) and it was noted that external shock synchrotron 
emission may be the strongest afterglow GeV component 



(see, e.g. [I^n et al |2008l Fan fc Piran|2008[ |Zou, Fan fc 



Piran|[2005 ) , the observation that this may be the dom- 



mant GeV emission over the whole burst, including the 

prompt phase, were surprising. 

Fo llowing t hese observations Ku mar fc Barniol Du-| 
(|2009ajb|) proposed a revolutionary model in which 



ran 



they revise a critical component of the standard external 
shock scenario. They suggest that there is no magnetic 
field amplification beyond the usual shock compression. 
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Namely, the downstream (shocked) magnetic field is just 
4r (where F is bulk Lorentz factor behind the shock) 
times the upstream circum burst field. In doing so they 
are able to fit the overall afterglow spectrum (ranging 
from optical to GeV), as the low magnetic field in the 
emitting region quenches the lower energy emission. Ad- 
ditionally they get rid of a n agging theoretic al problem - 
how are the fields amplified ( Gruzinov]|2001 )? 

The magnetic field plays a triple role in the 
synchrotron-shock acceleration mechanism. It acceler- 
ates the electrons and confines them to the shock, while 
they are accelerated and it also controls the synchrotron 
emission. A weaker magnetic fields poses two challenges 
to the model: cooling and confinement. First, a com- 
parison of the acceleration and the cooling times sets an 
absolute limit on the energy of synchrotron photon in 
the radiating fluid frame. Together with the hydrody- 
namics of the decelerating blast wave this puts a time 
dependent limit on the maximal energy of observed syn- 
chrotron photons. Photons above this limit are (almost 
certainly) not emitted by external shock synchrotron. Ef- 
flcient cooling poses another limit on the model. A sig- 
nificant (though not dominant) amount of energy is emit- 
ted in the GeV emission. This implies that the emitting 
electrons must be fast cooling. Otherwise the system 
would be inefficient and the energy requirement unrea- 
sonable. As the cooling takes place mostly in the down- 
stream region this last condition constrains the magnetic 
field there. While our original motivation was to examine 
the "unamplified" magnetic field scenario our analysis is 
more general and we allow for an amplification factor. 
We show that the observations of a significant GeV flux 



2 



Piran and Nakar 



poses strong limits on the downstream magnetic field. 
These limits can be translated to limits on the upstream 
circum burst field (in the case of no amplification) or on 
the amplification factor. 

Confinement is most important in the upstream re- 
gion, where the magnetic field is weakest. Thus, ob- 
servations of GeV photons limit the upstream magnetic 
field with a weak dependence on field amplification in the 
shock. Finally we turn to the influence of Inverse Comp- 
ton (IC) cooling on the observed GeV emission. Given 
the strong low energy (IR-optical) radiation fields (from 
the prompt, reverse shock and the forwards shock itself), 
the magnetic field density should be strong enough in or- 
der that IC cooling won't quench the GeV emission. This 
sets yet another, independent, limits on the downstream 
magnetic field. These considerations shed a direct light 
on the magnetic fields which are among the most elu- 
sive parameters of the external shock model. Note that 
here we derive constraints assuming that the external 
shock is adiabatic. If i t is radiative (as suggested e.g., by 
Ghisellini et al. 1120101 then the constraints will be more 



stringent. 

We examine in this letter these limits that arise from 
the GeV emission. We don't attempt to provide a com- 
plete solution to the whole multiwavelength afterglow. 
As such our analysis is very general and it depends only 
on the assumptions of synchrotron process and the blast 
wave hydrodynamics. 

2. THE MAXIMAL SYNCHROTRON ENERGY 

It is well known (e.g., dc Jage r fc: Harding|[l992 Lyu- 
tikov||2009[ |Kirk fc Reville,.2010] that by equatingTEe 
synchrotron cooling to the acceleration rate, one can ob- 
tain an upper bound on the maximal energy of a syn- 
chrotron photon: 



nipC 



Ma 



(1) 



This limit is in the fluid's rest frame, denoted by " ' ". 
h is Planck constant, nie is the electron rest mass, c is 
the speed of light and a is the flne structure constant. 
The maximal observed energy is larger by the Lorentz 
boost factor from the fluid frame to the observer frame, 
T) = [r(l — /3/i)(l -I- z)]~^, where F is the Lorentz factor 
of the fluid and fj. is cosine the angle between the fluid 
velocity and the line-of-sight. 

If the emission originates from a decelerating external 
shock then at any given time the observer receives si- 
multaneously photons emitted from a range of radii and 
therefore a range of Lorentz factors. The photons ob- 
served at a time t all satisfy t = T — Rfi/c, where T is the 
time since the explosion as measured in the source frame 
and R is the shock radius at time T. In a circum burst 
medium with mass density p (x this condition implies 
that along the line-of-sight (i.e., /i = 0) the Lorentz boost 
is Vios = 2r;o,/(l + z), where F^^^^ = i?,o,/4(4 - fc)ct. 
The Lorentz factor from any other direction, fj., (keeping 
t constant) is larger, but not necessarily the boost, which 
satisfles: 
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2(4 - fc) 



2(4-fc) 



(2) 



density medium(fc = 0), e.g., interstellar medium (ISM), 
T^max — 1-2 2?;os while in a circumburst wind (fc = 2) 
T^max = T^ios- For a self-similar adiabatic shock that 
propagates into ISM with a constant particle density, 
n, Tios,iSM = [llEfi/ie^TrmpC'nil + z)3]i/8 ^hile in 



a wind with a mass density p = Ar 



s.wind 



[9M/1287rc3A(l z)]i/4 dBlandford fc McKee| |l976 ) . 
Here E is the blast wave kinetic energy, rup is the proton 
mass and z is the burst's redshift. 

Using these expressions we flnd that the maximal ob- 
served energy of a synchrotron photon emitted by an 
adiabatic external shock during its decelerating phase is: 



9 GeV 
5 GeV 
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wind, 



(3) 



where we use the common notations of t2 — t/100, i?54 
£^/10^^ in c.g.s. units, etc. n is the ISM density and is 
the wind parameter in units of 5 • 10~^^ gr/cm~^. Note 
that the dependence on the burst's parameters of this 
rather general condition is very weak. Moreover, even 
when the burst redshift is unknown this limit peaks at 
z ~ 1. In principle this limit can be violated by syn- 
chrotron emission in a special magnetic fleld conflgura- 
tion, for example, if the electrons are accelerated by a 
weak magnetic field but radia te where the field is strong 
( |Lyutikov|2009l [Kumar l2010[ ) . Stfll, it is unlikely that af- 
tergiow photons 10 GeV are generated by synchrotron 
radiation. 

Such photons were already observed in two bursts. 
Most notable is GRB 940217, where an 18 GeV photon 
was observed 90 min after the burst by EGRET (iHurl^ 



et al. 



1994|) . An order of magnitude larger than the limTT 

ot Eq. [3|at this time, 2 GeV. Recently, LAT detected a 
33 GeV photon 82 s a fter the burst from GRB 090902B 



Maximizing T) for F > Vios shows that in a constant 



( Abdo et al. ||2009b[ ). The observation of one or two 
photons cannot rule out the external shock synchrotron 
model for GeV emission. Nevertheless, these observa- 
tions set a major difficulty to this model, especially if 
observations of late afterglow photons much above this 
limit will continue. 

3. CONFINEMENT AND COOLING 

If the observed GeV emission arises from synchrotron 
in an external shock it constrains both the upstream and 
the downstream magnetic fields. It constrains the first by 
the requirement that the radiating electrons are confined 
to the system. It also constrains the latter by the re- 
quirement that the radiating electrons radiate efliciently 
and therefore are cooling fast. The limits are particu- 
larly interesting in the case that the downstream fleld is 
not amplified beyond the usual shock compression and 
Bd = 4Fi3„. Such a case is interesting for this scenario 
because a rather weak magnetic field is needed in order 
to suppress the synchrotron MeV emission of the for- 
ward shock ( Kumar fc Barniol Duran||2009a ) . However, 
in order to consider the most general case, and since IC 
emission may suppress the MeV emission without affect- 
ing the GeV luminosity, we introduce an amplification 
factor, /b > 1, so that Bd — AJbTBu- In the follow- 
ing we presents these constraints using Jb and upstream 
field, which for abbreviation we denote simply as B. 
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3.1. Confinement 

The accelerated electrons have to be confined to the 
shock region, otherwise they escape. This sets a second 
limit on the maximal synchrotron energy. The confine- 



fdR/r 



§d =4fBr§u 




Fig. 1. — The motion of electrons in the upstream (right - red) 
and downstream (left - purple). Note that the shock moves for- 
wards and slows down while the electron gyrates in the upstream 
magnetic field. 

mcnt criterion on the Gyro radius in the upstream of an 
electron with Lorentz factor 7^ is: 



Rl = 



Tj[jrie 
eB 



<fuTR 



(4) 



The external shock slows down while it propagates and 
it won't catch the rapid electron by the time it turned 
i?L/r but somewhat later (see Fig. [T]). This leads to 
the confinement factor, /u < 1. Detailed calculations 
yield ~ 1/3 for an ISM and 1/2 for a wind, we denote 
fu,x — fu/x. This last inequality leads to a maximal 
confined electron's Lorentz factor (in the shock's frame): 



^conf 



eB 



-JuR 



(5) 



For confinement in the downstream the Larmour radius 
should be compared with the shock's thickness fdR/^ 
where fd ~ 0.1. The resulting maximal electron Lorentz 
factor is Afdfs / fu times larger than the one given by Eq. 
[5] for the upstream. Since the electron has to be confined 
in both upstream and downstream the comparable or 
smaller former limit is the critical one. 

The synchrotron frequency of electrons in the down- 
stream is larger by a factor of A/b than in the upstream. 
Therefore the confinement upper limit for synchrotron 
external shock photons is: 



hv, 



conf 



eBY flR^T^ 



mc J (1 + z)c2 



(6) 



The average Lorentz factor and the blast wave radius 
in a constant density environment that correspond s to 
an observer time t are ( Sari, Piran fc Narayan|[l998 ) : 

1/8 /i + ,x3/8 



180 



-3/8 



(7) 



and 



i? = 2 • 10^^ cm 



1/4 



1 I \ -1/4 

'^'^ tT- (8) 



n ) V 2 

Substitution of these expressions into Eqs. [5]and|6]yields: 
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3/4 



1 I \ -3/4 
1 + Z\ 1/4 

'■2 



(10) 

These equations imply that if the observed GeV photons 
were generated by an external shock synchrotron the up- 
stream magnetic field must satisfy: 



B>20fiG 



10 GeV 



Turning now to a wind environment the average 
Lorentz factor and blast wave radius that cor responds 
to an observer time t are ( |Chevaher fc"Li|[2000l ) : 



£^54 



1/4 



X 1/4 
Z \ ,-1/4 



(12) 



and 

R^u^nd = 3.7 X lO^^ Cm 



E54 

~a7 



1/2 



1 I \ -1/2 
1 + 2 \ 1/2 

'■2 



A substitution of these expression into Eqs. and mi 
yields: 



I'conf-wind — 10^/u,i-^-5 ^^4~^ 2 



hv, 



con f —wind 
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l + z 



X -1/2 
Z \ 1/2 

'■2 

(14) 

_ 3 

r 

(15) 

The lower Loretnz factor and the smaller radius (at early 
time) as compared to an ISM leads to a more stringent 
constrain on the upstream magnetic field of a wind: 



E^ 



B>70hG 



10 GeV 



E^i 
n 



l + z 



(16) 



3.2. Cooling 

The conditions found in the previous subsection are 
sufficient to produce the highest energy photons. How- 
ever, the bulk of the energy typically observed in the LAT 
comes at lower energies, around 100 MeV. To obtain an 
efficient, fast cooling, 100 MeV emission one needs the 
cooling frequency, Vc < 100 MeV. In both the upstream 
and downstream frames the maximal energy, VMax is ob- 
tained from the condition tcooi = tacc, while v^onf is 
obtained from the condition tdyn = tacc- Since cooling 
is more efficient in the downstream, comparing the two 
conditions in the downstream frame implies that Vc (for 
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which tr 



tdyn) satisfies 



'Ma 



^conf,d 



(17) 



where i^confM = {^IdfB/fuYvconf- This equation of the 
cooling frequency is, of course, equivalent to the one ob- 
tained by the traditional derivation (e.g.. Sari, Piran & 
Narayan 1998). Thus, instead of calculating the exact 
value 01 fd we use t he normalization of derived by 



Granot k. Sari (2002) finding 



60 GeV {fBB-^)-^"'' 
200 GeV (/bS„5)-' 



to 



g54 

A, 



2 ' ISM 
wind 

The wind constrain is again stronger since v^onf is lower 
while VMax is independent of the environment. The con- 
dition hvc < lOOMeV, will impose now: 



fBB> 



85AiG (^)"'/'(i±^)'/'t2-^/^ ISM 



125 ^G {^] (i±^f /^r^/^^ 



A, 



wind 



(19) 

This is a relatively high upstream field and it suggests 
that either the external density is extremely low or that 
there is at least some magnetic field amplification at the 
shock. 

The system has some interesting properties when there 
is no magnetic field amplification. In such a case the 
upstream and downstream timescales are comparable to 
within an order of magnitude and each electron is char- 
acterized by only two time scales, its acceleration time, 
tacc and its cooling time, tcooi- It is the interplay be- 
tween these two time scales and the dynamical time scale, 
tdyn J together with the requirement that the electron is 
bound to the system, that determine the critical electron 
Lorentz factors, jMax, Iconf, and 7c. The first satisfies 
taccilMax) = tcooi{lMax)- It is independent of the sys- 
tem dynamical time and as it turns out VMax is also 
independent of the strength of the magnetic field (Eq. 
fll). Next, equating tacc = tdyn results in the maximal 
Torentz factor an electron can achieve, when cooling and 
escape are ignored. In our case, an electron that satis- 
fies tacc — tdyn is also the maximal energy of an electron 
that is confined to the system, thus, tacdlconf) = tdyn- 
Finally, 7c is the Lorentz factor of the electron that cool 
over the dynamical time, i.e., tcooiilc) = tdyn- Therefore 
there are two possibihties. If tacdlMax) = tcooi{lMax) < 
tdyn then the maximal electron's energy is limited by 
cooling, which takes place on shorter time than tdyn- 
Thus, less energetic electrons than jMax can cool over 
tdyn implying i/c < VMax < i^conf- In such case a 
break in th espectrum is observed at Vc and a cutoff at 
I'M, while Vconf is not observed. On the other hand if 
tdyn < taccilMax) = tcooiijMax) then the maximal elec- 
tron's energy is limited by the ability to accelerate and 
confine the electrons over tdym where cooling do not play 
any role. Thus, i^conf < VMax < Vc where ^Max and 7c 
electrons do not exist in the system (cannot be acceler- 
ated) . Therefore a cutoff is observed at Vconf and Vc and 
VMax are not observed. For an ISM there is a triple coin- 
cidence of the three energies for our canonical parameters 
when B = 20/LtG so i^conf = VMax = Vc- 



4. INVERSE COMPTON 

So far we have ignored the possible effects of Inverse 
Compton (IC) cooling on the synchrotron emitting elec- 
trons. However, the strong lower energy radiation fields 
may lead to IC cooling whose effect we consider now. 
This radiation field may be the external forward shock 
synchrotron field, but it may also be any other source. 
The IC cooling may reduce the maximal energy bellow 
maximal synchrotron energy (Eq. [T]). It leads to an im- 
plicit equation for the new (lower) maximal energy: 



(20) 



where Y is the Compton parameter in the downstream 
(where the dominant cooling takes place) and is 
the Lorentz factor of the electrons whose synchrotron 
frequency is v,^, j^- Clearly, for efficient GeV emission 
^{1m ic) must be smaller than unity and in this case 
1m IC ~ iMax and we recover Eq. [ij In the following we 
assume that F < 1 and ^^.ic ~ iMax and consider what 
are the conditions for this to hold. 

As "/Max is extremely high IC scattering is mostly in 
the Klein Nishina (KN) regime and is ineffective. Still, 
the radiation field (in the shock frame), U^ad' below the 
relevant KN frequency, j^kjv(7)' 



hvj,^ (7) 



7(1 + ^) 



(21) 



may be large enough leading to an effective Compton 
parameter: 



Y{-/Max) 



,{lMax)) 



(22) 



One can estimate U'^a^ using a given model for the ex- 
ternal shock emission, but given the multi wavelength 
observations at early time it is much more useful to con- 
strain it directly from the observations. This is especially 
important since, at least in the earlier part of the after- 
glow, other radiation fields (e.g. prompt emission and 
reverse shock emission) co-exist with the external shock 
and this is the best way to incorporate their contribu- 
tions. Given an observed spectral flux density 



U'rad{v < ''kn) = max{z^F^(iy < i^^.„)} 



t(l -I- z) 



i?3 ' 

(23) 

where d is the comoving distance. In the follow- 
ing we measure F^, in units of mJy (i.e., F^, — 
i^26lO~^^ergs/cm^Hz). ma.x{i'F^{i' < i'k^)} is the max- 
imal value of vFi, at frequencies below . Below we 
assume that this maximum is at v^^^^ , an assumption that 
can be tested, and cor rect ed, given the observations. 

[Tslwe find that v^^ril Max) is 



Using Eqs. j7)[|)[T|and 

typically in the IK when t 
cation: 



icre is no strong field amplifi- 



r ilMax) 



0.05eVv/7^Sr;(f4)^« (i±2) 1^ i2 - ISM, 
0.02 eVy/s^ 



'1 + 2^ 

V 2 J 



wind. 



(24) 
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The relevant Y parameter assuming i'F^{< peaks 



at 



y{lMax) = 



0.2 F26d28'(/BB-5) 
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Additionally, in order to produce efficiently the emission 
observed by LAT at > 100 MeV, the corresponding elec- 
trons cannot be cooled effectively via IC, and the Comp- 
ton parameter should be small. The relevant KN fre- 
quency is: 

(100 MeV) = 
0.5eV (/si?-5)^ (^)^^(i^^)"^i2-J ISM, 

0.2 eV (/bB_5)5 (f^) ' i^y' wind. 

(26) 

and assuming vF^{< i^^^) peaks at v^^^^ the Compton 
Y is: 



y(100 MeV) 



1.3F26d28'(/i3S_5)-i 
180 F26d28'(/BB-5)-^ (ff) 



4,-1-4, 



2 1+z 



t| ISM, 



2 ^2 



wind. 



(27) 

Therefore, if there is no strong field amplification a mod- 
est IR-Optical flux of a few mJy for an ISM and 10 /iJy 
for a wind is enough to suppress the GeV flux. There- 
fore a clear test of the external shock synchrotron GeV 
emission, without field amplification is that we should 
not observe a strong optical flux at the time when we 
observe GeV emission. A simultaneous detection of a 
strong optical emission and a GeV emission will rule out 
this model. 

Although we do not consider here a complete model for 
the external shock synchrotron emission it is important 
to note that when such model is constructed KN effects 
should be carefully considered. The reason is that if the 
magnetic field is not amplified close to equipartition level 
then the radiation energy density is much larger than the 
magnetic field energy density in the downstream frame. 
IC by electrons emitting the MeV-GeV emission will cer- 
tainly be in the KN regime, but over a large range of the 
parameter phase the Y parameter of these fast cooling 
electrons (which now depends on the electron energy) will 
be larger than 1. In which ca se the synchrotron spectrum 
will be altered significantly (Nakar, Ando & Sari 2009). 



An extreme example to a KN effect may be apparent if 
there is no field amplification at all. In such case j^^^^ 
of the fast cooling electrons is below the the typical syn- 
j^^tron frequency VMax-^i their Compton parameter is 
smaller than unity and the synchrotron emission is not 
■vetedugly affected. But, once Y becomes larger than 1 for 
fast cooling electrons, the back-reaction of the IC emis- 
sion on the electrons distribution results in the cooling 
frequency "jumping" on a short time scale by orders of 
magnitude, signific antly revising the whole synchrotron 
spectrum ( |Nakar, Ando fc Sari||2009[ ) . 



5. CONCLUSIONS 

Using the confinement and cooling conditions we have 
obtained limits on the values of the magnetic fields 
needed in the downstream and upstream regions in order 
to produce the observed GRB GeV emission via an ex- 
ternal shock synchrotron. These constrains are based on 
minimal assumptions of synchrotron cooling and blast 
wave hydrodynamics. Both are essential ingredients of 
the external shock synchrotron model. The arguments 
we present allow us to explore directly the magnetic fields 
in both upstream and downstream regions, which are 
among the least constrained physical parameters of the 
model. 

We find that with no amplification the minimal fields 
required are on the high end (lOO/i G), unless the external 
density is very low. The limits are even higher for a 
radiative solution. It is, of course, possible that this is 
a condition for GeV emission. However, the detection 
of GeV emission from all MeV bright GRBs that are 
within the LAT viewing angle suggests that the emission 
is generic. In this case at least a modest amplification is 
probably needed. 

Finally, we point out two critical predictions of the ex- 
ternal shock synchrotron model: (i) No detection of late 
very energetic (> 10 GeV) photons and (ii) No simultane- 
ous detection of a bright (>mJy) IR-optical (depending 
on the specific case) signal with the GeV photons un- 
less the upstream magnetic field is strongly amplified in 
the shock. Continued observations should be compared 
with these predictions and can provide future tests of 
this model. 
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